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Abstract 

Liquids in contact with solids are submitted to intermolecular forces inferring den- 
sity gradients at the walls. The van der Waals forces make liquid heterogeneous, the 
stress tensor is not any more spherical as in homogeneous bulks and it is possible to 
obtain stable thin liquid films wetting vertical walls up to altitudes that incompress- 
ible fluid models are not forecasting. Application to micro tubes of xylem enables 
to understand why the ascent of sap is possible for very high trees like sequoias or 
giant eucalyptus. 

Key words: nanofilms, inhomogeneous liquids, van der Waals forces, ascent of sap, 
high trees. 

PACS: 68.65.k, 82.45.Mp, 87.10.-he, 87.15.Kg, 87.15.La 



1 Introduction 

In Amazing numbers in biology, Flindt reports an eucalyptus of 128 meters 
and a giant sequoia of 135 meters pj. However, biophysical determination of 
maximum size to which trees can grow is not well understood [2]. A main 
problem with the understanding of tall trees is why the sap is able to reach 
so high levels. 
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Xylem tube diameters range between 50 and 400 /im; the crude sap contains 
diluted salts but its physical properties are roughly those of water. Conse- 
quently, hydro dynamic, capillarity and osmotic pressure create a sap ascent of 
only few tens of meters To explain the sap ascent phenomenon, Dixon and 
Joly proposed in 1894 a cohesion-tension model [1], followed by a quantitative 
attempt by van der Honert in 1948 [5]: liquids may be subjected to tensions 
generating negative pressures compensating gravity effects. Nevertheless, ther- 
modynamic states are strongly metastable and can generate cavitation causing 
embolisms in xylem tubes made of dead cells 

As pointed out in Ref. [7], a turning-point in the pro and con debate on the 
sap ascent was the experiment of Preston in 1952 who demonstrated that tall 
trees survived overlapping double saw-cuts made through the cross-sectional 
area of the trunk to sever all xylem elements [8J. This result, confirmed later 
by several authors (e.g. Mackey & Weatherley in 1973 [9J; Eisenhut, in 1988 
[TU] : Benkert et al in 1991 [H]), was obviously not in agreement with the 
cohesion-tension theory. Using a xylem pressure probe. Balling and Zimmer- 
mann showed up that, in many circumstances, this apparatus does not measure 
any water tension [12] . Since these experiments, Zimmermann et al questioned 
the cohesion-tension theory in Ref. [13]: xylem tension exceeding 0.6 Mpa was 
not observed and in normal state most vessels were found to be embolized at 
a level corresponding to sixty meter high |14p1: consequently trees growing 
taller than a few tens of meter range are not foreseeable. Recently in their re- 
view article Zimmermann et al demonstrate that the present interpretation of 
the pressure bomb data is based on a misconception and that negative xylem 
pressure values of several megapascals do not exist since xylem sap composi- 
tion, the features of the xylem wall and the hydraulic coupling of the xylem 
with the tissue prevent the development of stable tensions larger than about 
1 MPa. Moreover, gas-vapor transportation in xylem tubes was found at the 
top of height trees [7] . 

In this paper, our ambition is to present an understanding of the ascent of 
sap in very high trees, different from the cohesion-tension theory: at a higher 



It is interesting to remark that xylem tube diameters range between 50 and 400 
fim and the crude sap is a liquid bulk with a superficial tension a lower than the 
superficial tension of pure water which is 72.5 cgs at 20° Celsius. Let us consider 
a microscopic gaz- vapor bubble inside the crude sap with a diameter 2R smaller 
than xylem tube diameters. The difference between the gaz-vapor pressure and the 
liquid sap pressure can be expressed by the Laplace formula: — Pi = 2a/ R. 
The vapor-gas pressure is positive and omitted with respect to \Pi\; consequently 
unstable bubbles appear when R > —2a/Pi . For a negative pressure Pi = —0.6 
MPa in the sap as pointed out by experiments, we obtain R > 0.24 //m. Due to 
the diameter range of xylem tubes and the diameter range of bubbles, the Laplace 
formula is valid at equilibrium for such bubble sizes [15J; dynamical bubbles appear 
spontaneously from germs naturally included in the liquid when the tubes are filled 
with the crude sap and cavitation makes the tubes embolized. 
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level than a few tens of meters - corresponding to the pulling of water by 
hydrodynamic, capillary and osmotic pressure - we assume that tubes may 
be embolized. In addition, we assume also that a thin liquid film - with a 
thickness of a few nanometers [TB|17lll8j - wets xylem walls up to the top of 
the tree. At this scale, long range molecular forces stratify liquids [19] and the 
ratio between tube diameter and sap film thickness allows us to consider tube 
walls as plane surfaces; consequently the problem of sap ascent in vertical 
tubes is similar to the rise of a liquid film damping a vertical plane wall. 
The sap motion in xylem tube can be suitably explained by the transpiration 
through micropores located in tree leaves PI7] : evaporation changes the liquid 
layer thickness implying driving of sap as explained in Ref. [20] . 
Consequently, this paper aims to prove, in non-evaporating case, the existence 
at equilibrium of thin films of water wetting vertical walls up to a same order 
of altitude than the height of very tall trees. 

The recent development of experimental technics allows us to observe physi- 
cal phenomena at length scales of a few nanometers. This nanophysics reveals 
behaviors often surprising and basically different from those that can be ob- 
served at a microscopic scale [21j . 

At the end of the nineteenth century, the fluid inhomogeneity in liquid-vapor 
interfaces was taken into account by considering a volume energy depending 
on space density derivative [22|23] . In the first part of the twentieth century, 
Rocard obtained a thermodynamical justification of the model by an origi- 
nal step in kinetic theory of gases [24j. This van der Waals square-gradient 
functional is unable to model repulsive force contributions and misses the 
dominant damped oscillatory packing structure of liquid layers near a sub- 
strate wall [25)126] . Furthermore, the decay lengths are only correct close to 
the liquid-vapor critical point where the damped oscillatory structure is sub- 
dominant [27P5] . Recently, in mean field theory, weighted density-functional 
has been used to explicitly demonstrate the dominance of this structural con- 
tribution in van der Waals thin films and to take into account long-wavelength 
capillary-wave fluctuations as in papers that renormalize the square-gradient 
functional to include capillary wave fluctuations [29|30] . In contrast, fluctua- 
tions strongly damp oscillatory structure and it is mainly for this reason that 
van der Waals' original prediction of a hyperbolic tangent is so close to simu- 
lations and experiments [21]. 

To get an analytic expression in density-functional theory for a thin liquid 
film near a solid wall, we add a liquid density-functional at the solid surface 
to the square-gradient functional representing closely hquid-vapor interface 
free energy. This kind of functional is well-known in the literature, as the gen- 
eral background studied by Nakanishi and Fisher [32] ■ It was used by Cahn 
in a phenomenological form, in a well-known paper studying wetting near a 
critical point [33]. An asymptotic expression is obtained in \3^ with an ap- 
proximation of hard sphere molecules and London potentials for liquid-liquid 
and solid-liquid interactions: by using London or Lennard- Jones potentials, we 
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took into account the power-law behavior which is dominant in a thin hquid 
film in contact with a solid. In this paper, the effects of the vapor bulk bor- 
dering the liquid film are simply expressed with an other density-functional of 
energy located on a mathematical surface as a dividing-like surface for liquid- 
vapor interfaces of a few Angstrom thickness. 

With this functional, we obtain the equations of equilibrium [35] and bound- 
ary conditions [36] for a thin vertical liquid film damping a vertical solid wall 
and we can compute the liquid layer thickness as a function of the film level. 
Moreover, the normal stress vector acting on the wall is constant through the 
liquid layer and corresponds to the gas-vapor bulk pressure which is currently 
the atmospheric pressure; no negative pressure appears in the liquid layer. 
As in [37], several methods can be used to study the stability of a thin liq- 
uid film in equilibrium. In our case, the disjoining pressure of a thin liquid 
layers is a well adapted tool for very thin films. By using Gibbs free energy 
per unit area for the liquid layer as a function of the thickness, we are able 
to obtain the minimal thickness for which a stable wetting film damps a solid 
wall. The minimal thickness is associated with the pancake layer when the 
film is bordering the dry solid wall [T6|38|39j and corresponds to a maximal 
altitude. Numerical calculations associated with physical values for water yield 
the maximal film altitude for a silicon wall and a less hydrophile material. In 
the two cases, we obtain an approximative maximum level corresponding to a 
good height order for the tallest trees. 



2 Definition and well-known results of the disjoining pressure 

Without redoing or demonstrating the main results of Derjaguin et al [TB] 
related to thin liquid films and the well-known disjoining pressure, we enu- 
merate the properties we apply in the problem of rise of a liquid on a vertical 
wall. 

In this paper, we consider fluids and solids at a given temperature 9. The film 
is thin enough such that the gravity effect is neglected across the liquid layer. 
The hydrostatic pressure in a thin liquid layer included between a solid wall 
and a vapor bulk differs from the pressure in the contiguous liquid phase. At 
equilibrium, the additional interlayer pressure is called the disjoining pressure. 
Clearly, a disjoining pressure could be measured by applying an external pres- 
sure to keep the complete layer in equilibrium. The measure of a disjoining 
pressure is either the additional pressure on the surface or the drop in the 
pressure within the mother hulks which produce the layer. In both cases, the 
forces arising during the thinning of a film of uniform thickness h produce the 
disjoining pressure n(/i) of the liquid layer with the surrounding phases; the 
disjoining pressure is equal to the difference between the pressure P^^ on the 
interfacial surface (which is the pressure of the vapor mother hulk of density 
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p^J and the pressure Pb in the liquid mother bulk (density p^) from which 
the hquid layer extends (this is the reason for which Derjaguin used the term 
mother bulk [TB], page 32) : 



U{h) 



The most classical apparatus to measure the disjoining pressure is due to 
Sheludko [10] and is described on Fig. (1). Let us consider the Gibbs free 
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Fig. 1. Diagram of the technique for determining the disjoining pressure isotherms 
of wetting films on a solid substrate: a circular wetting film is formed on a flat 
substrate to which a microporous filter is clamped. A pipe connects the filter filled 
with the liquid to a reservoir containing the liquid mother hulk that can be moved 
by a micrometric device. As we will see in section 5, the thickness h of the film 
depends on H in a convenient domain of H values where the wetting film is stable. 
The disjoining pressure is equal to 11 = (p^ — p^^)gH, where g is the acceleration 
of gravity (From Ref. JTW . page 332). 

energy of the liquid layer (thermodynamic potential). As proved by Derjaguin 
et al in Ref. (p^, Chapter 2), the Gibbs free energy per unit area G can be 
expressed as a function of h : 



dG{h) 
dh 



Eq. ([2]) can be integrated as 



G{h) 



-Yi{h). 



n(/i) dh 



(3) 



where /i = is associated with the dry wall in contact with the vapor bulk 
and h = +oo is associated with a wall in contact with liquid bulk when the 
value of G is 0. 

An important property related to the problem of wetting is associated with 
the well-known spreading coefficient : 



S 



7c 



7. 
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where Igy^lg^yli^y are respectively the solid- vapor, solid-liquid and liquid- 
vapor free energies per unit area of interfaces. The energy of the liquid layer 
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per unit area can be written as 



When h = 0,we obtain the energy 7^^ of the dry sohd wall; when h 
obtain 7^^ +7^^- complete wetting of a liquid on a solid wall, the spreading 
coefficient is positive and the Gibbs free energy G looks like in Fig. 2. 



Fig. 2. The construction of the tangent to curve G(h) from point W of coordinates 
(0,G(0)) involves point P; point W is associated with a high-energy surface of the 
dry wall and point P is associated with the pancake thickness hp where the film 
coexists with the dry wall; it is the smallest stable thickness for the liquid layer. 

The conditions of stability of a thin liquid layer essentially depend on phases 
between which the film is sandwiched. In case of a single film in equilibrium 
with the vapor and a solid substrate, the stability condition is : 

The coexistence of two film segments with different thicknesses is a phe- 
nomenon which can be interpreted with the equality of chemical potential 
and superficial tension of the two films. A spectacular case corresponds to 
the coexistence of a liquid film of thickness hp and the dry solid wall associ- 
ated with /i = 0. The film is the so-called pancake layer corresponding to the 
condition : 

G{{))=G{hp) + hpIl{hp). (5) 
Eq. ([5]) expresses that the value of the Legendre transformation of G{h) at hp 
is equal to G'(O). Liquid films of thickness h > hp are stable and liquid films 
of thickness h < hp are metastable or unstable. 



3 The study of inhomogeneous fluids by using a square- gradient 
approximation and surface-density functionals at bordering walls 



The modern understanding of liquid-vapor interfaces begins with papers of 
van der Waals and the square-gradient approximation for the free energy of 
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inhomogeneous fluids. In current approaches, it is possible to give formal exact 
expressions of the free energy in terms of pair-distribution function and the 
direct correlation function [H]. In practice, these complex expressions must 
be approximated to lead to a compromise between accuracy and simplicity. 
When we are confronted with such complications, the primitive mean-field 
models are generally inadequate and the obtained qualitative picture is no 
more sufficient. The main alternatives are density-functional theories which 
are a lot simpler than the Ornstein-Zernike equation in statistical mechanics 
since the local density is a functional at each point of the fluid [TTf^ . We 
use this rough approximation enabling us to compute analytically the density 
profiles of simple fiuids. Nevertheless, we take into account surface effects 
and repulsive forces by adding density-functionals at boundary surfaces. The 
density-functional of the inhomogeneous fluid in a domain O of boundary dO 
is taken in the form : 



e dv + w ds. (6) 

o J JdO 

The flrst integral is associated with a square-gradient approximation when we 
introduce a speciflc free energy of the fluid at a given temperature 9, 

e = e{p,(3) 

as a function of density p and (3 = (grad p)^. Speciflc free energy e characterizes 
together the fluid properties of compressibility and molecular capillarity of 
liquid-vapor interfaces. In accordance with gas kinetic theory, 

= pa(p) + ^ (grad p)^ (7) 

where term (A/2) (grad p)^ is added to the volume free energy pa{p) of a 
compressible fluid and A = 2pe'p{p,(3) is assumed to be constant at a given 
temperature [21] . Speciflc free energy a enables to connect continuously liquid 
and vapor bulks and pressure -P(p) = p^a'^^p) is similar to van der Waals one. 

Near a solid wall, London potentials of liquid-liquid and liquid-solid interac- 
tions are : 

'^11 = — ^ ! when r > ai and = oo when r < ai , 
V'ls — ~~j ; when r > 6 and (pi^ = oo when r < 6 , 

where cu and cis are two positive constants associated with Hamaker constants, 
0"; and as denote fluid and solid molecular diameters, 5 = |( cr;+ cTs) is the 
minimal distance between centers of fluid and solid molecules. Forces between 
hquid and solid are short range and can be described simply by adding a spe- 
cial energy at the surface. This energy is the contribution to the solid/fluid 
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interfacial energy which comes from direct contact. This is not the entire inter- 
facial energy: another contribution comes from the distortions in the density 
profile near the wall [S51l34f38j . For a plane solid wall (at a molecular scale), 
this surface free energy is in the form : 

0(P) = -7iP+^72P'- (8) 

Here p denotes the fluid density value at surface (5); constants 7^^, 72 and A 
are positive and given by the mean field approximation : 

where mi and denote respectively masses of fluid and solid molecules, p^^i 
is the solid density [51] . 

We consider a plane liquid layer contiguous to its vapor bulk and in contact 
with a plane solid wall (5); the z-axis is perpendicular to the solid surface. 
The liquid film thickness is denoted by /i; the conditions in the vapor bulk are 
gradp = and Ap = with A denoting the Laplace operator. 
Far below from the critical point of the fluid, a way to compute the total free 
energy of the complete liquid-vapor layer is to add the surface energy of the 
solid wall (S) at z = 0, the energy of the liquid layer (L) located between z = 
and z = h, the energy of the sharp liquid-vapor interface of a few Angstrom 
thickness assimilated to a surface (S) a.t z = h and the energy of the vapor 
layer located between z = h and z = +00 [12]. The liquid at level z = h is 
situated at a distance order of two molecular diameters from the vapor bulk 
and the vapor has a negligible density with respect to the liquid density [13]. 
In our model, the two last energies can be expressed with writing a unique 
energy ip per unit surface located on the mathematical surface (S) at z = h : 
by a calculus like in Ref. [M], we can write ip in the same form than Rel. ([8]) 
and also expressed as in Ref. [38] in the form ipi^p) = — 75P + 1 74 p^; but with 
a wall corresponding to a negligible density, 75 ~ 0, the surface free energy ip 
is reduced to : 

^(P) = y P^ (10) 

where p is the liquid density at level z = h and 74 is associated with a dis- 
tance d of the order of the fluid molecular diameter (then (i ~ 5 and 74 ~ 72). 
Consequently, due to the small vapor density, the surface free energy ip is the 
same than the surface free energy of a liquid in contact with a vacuum. 
Complementary to this argumentation, we will see, in section 4, that the 
boundary condition at surface (S) associated with surface energy (ITO!) yields 
a density value corresponding to an intermediate density between liquid and 
vapor and which can be considered as a density value of a dividing-like surface 
separating liquid and vapor inside the liquid-vapor interface. 



8 



Density-functional ([6]) of the hquid-vapor layer gets the final form : 
F = [ [ [ e dv + [ [ (p ds + [ [ ip ds. 

J J J(L) J J{S) J J(S) 



4 Equation of equilibrium and boundary conditions of a thin liquid 
layer contiguous to its vapor bulk and in contact with a vertical 
plane solid wall 

In case of equilibrium, functional ffTTl) is stationary and yields the equation of 
equilibrium and the boundary conditions [5B|37ll43j . 



4-1 Equation of equilibrium 



The equation of equilibrium is : 

div (T — p grad Q = 0, 
where Q is the body force potential and cr the stress tensor generalization 

mag, 

(T = — pi — X grad p ® grad p, 
with p = p^e'p — p div (A grad p). 

Let us consider an isothermal vertical film of liquid bounded respectively by 
a flat solid wall and a vapor bulk; then 



div (T + pg\ = Q 



:i2) 



in orthogonal system, where i is the downward direction of coordinate x (the 
gravity potential is f2 = —gx). 

The coordinate z being external and normal to the flat vertical solid wall, 
spatial density derivatives are negligible in directions other than direction of z 
corresponding to a very strong gradient of density of the liquid normally to the 
layer and a weak inhomogeneity along the film. In the complete liquid-vapor 
layer (we call interlayer) , 
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and Eq. (fT2|) yields a constant value at level x for the eigenvalue 03, 



where P^.^^ denotes the pressure P{Py^ ) in the vapor bulk of density p^^ 
bounding the liquid layer at level x. In the interlayer, eigenvalues 01,02 are 
not constant but depend on the distance z to the solid wall. In all the fluid, 
Eq. f[T2|) can also be written [35] : 

grad ( /i — A Ap — (yf X ) = , (13) 

where /i is the chemical potential (at a temperature 6) defined to an unknown 
additive constant. We note that Eqs. fll2til3p are independent of surface ener- 
gies P and (fTO|). 

The chemical potential is a function of P (and 6); due to the equation of state 

for pressure P, the chemical potential can be also expressed as a function of 

p (and 9). We choose as reference chemical potential p^ = Po{p) null for bulks 

of densities pi and p^ of phase equilibrium. Due to Maxwell rule, the volume 

free energy associated with p^ is go{p) — Po where Po = P{pi) = P{Pv) is the 

fp 

bulk pressure and go{p) = / Poip) dp is null for the liquid and vapor bulks 

Pv 

of phase equilibrium. The pressure P is : 

P{p) = pp,{p)-go{p) +Po. (14) 
Thanks to Eq. f|T3l) . we obtain in all the fluid and not only in the interlayer : 

Poip) - ^^P -9 x = Poipb), 

where Po{Pb) is the chemical potential value of a liquid mother bulk of density 
Pfj such that PoiPb) = PoiPvJ^ where p^^ is the density of the vapor mother 
bulk bounding the layer at level x = 0. This property is due to Eq. f|T3|) which 
is valid not only in the liquid but also in all the fluid independently of the 
surface energies in the density-functional flTT]) . Equation f[T^ is also valid in 
the sharp liquid-vapor interface. 

We must emphasize that P{pi) and P{p^^) are unequal as for drop or bubble 
bulk pressures. Likewise, we define a liquid mother bulk of density p^^ at level 
X such that Po{,PbJ = PoiPv^J with PipbJ ^iPv^J'^ Pb^ is not a fluid density 
in the liquid layer but density in the liquid bulk from which the interlayer can 
extend. Then, 

^ ^P = Poip) - PoiPbJ with p^{pf,J=p^{pf,)+gx (15) 
and in the interlayer 



d'^p 

A ^ = PbAP)' with /ife^(p) = p^{p) - p^iPhJ (16) 
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4-2 Boundary conditions 



The condition at the sohd wall (5) associated with the free surface energy ([8]) 
yields [36] : 

where n is the external normal direction to the fluid. Eq. (fT7I) yields : 




The condition at the liquid- vapor interface (S) associated with the free surface 
energy (ITU]) yields : 



\z = h 



As we will see in Sects. 5 and 6, Rel. ffTS]) takes into account the density 
aX z = h which is smaller, but of the same order, than liquid density. Due 

to the numerical values of A and 74 in Sect. 6, the density derivative 

dz 

is large with respect to the variations of the density in the interlayer and 
corresponds to the drop of density in the liquid-vapor interface in continuous 
model. Consequently, Rel. f|T8l) defines the film thickness by introducing a 
reference point inside the liquid-vapor interface bordering the liquid layer with 
a convenient density at surface z = h considered as a kind of dividing-like 
surface in a continuous model for the liquid- vapor interface (|31j, Chapter 3). 



5 The disjoining pressure for vertical liquid films 



Eq. ([T]) can be extended with the disjoining pressure at level x [16] : 

n = p,,^ - , 

where Pb^ and P„^^ are the pressures in liquid and vapor mother bulks cor- 
responding to level X. At a given temperature ^, 11 is a function of p^^ or 
equivalently a function of x. Let us denote by 

9bAp) = 9o{p) - 9o{pbJ - f^oiPbJip - PbJ, (19) 
the primitive of p^SP) ^^^^ Pb^- Consequently, from Eq. f[T^ . 

ii{pJ = -gbAp.J. (20) 
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and an integration of Eq. (fT6|) yields : 

= ^7,,(p) + n(p. ). (21) 



A dp' 
2 \di 



The reference chemical potential linearized near p, is Poip) — ~iP~Pi) where 

Pi 

ci is the isothermal sound velocity in liquid bulk pi at temperature 6 |l5j. In 
the liquid part of the liquid- vapor film, Eq. (|T6l) of density profile yields : 

>^-r^ = —{p - Pb) - 9 X = —{p - PbJ with pf,^ = pi, + ^ g X. (22) 
az pi pi Ci 

The reference chemical potential linearized near p^ is Po{p) = ~(P ~ Pv) 

Pv 

where is the isothermal sound velocity in vapor bulk p^ at temperature 9 
[15] . In the vapor part of the liquid- vapor film, 

^J^ = f^P-pJ-9^=f^P-PvJ with p^^^= p,^ + ^gx. 

p^ Py 

Due to Eq. ( |T5|) . Po{p) has the same value for p^^ and p^^, then 
— (Pb. - Pi) = PoiPbJ = PoiPv.) = —(Pv. - PJ> and 

Pi Pv 

p =p (i + i^p^^-pA. 

In the liquid and vapor parts of the interlayer we have, 

9o{p) = 7^{p - Pif (liquid) and go{p) = 7r^{p - P^f (vapor). 

'^Pl ^Pv 

From Eqs (fT9l) - (l20l) we deduce immediately the disjoining pressure at level x 

(23) 



^^P. + P.-pj2+^^(^--^')' 



n(P6j = ^(A-P6. 
^Pi 

Due to p, (2 + % Sp^^~P^^ \ ^ + p we get 
\ Pi / 



Pi 



n(p.J^|^(p?-pL). 



At level X = 0, the liquid mother bulk density is closely equal to p^ (density 
of liquid in phase equilibrium) and because of Eq. (122]) . 11 can be considered 
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as a function of x : 



n(a;) = -pigxil + j^ 



(24) 



Now, we consider a film of thickness hx at level x; the density profile in the 
liquid part of the liquid-vapor film is solution of : 



. , .dp , ^dp 

with A— = -7i + 72 p\^^^ and A— 

|z=0 W^l^^ft^ 



Quantities r and d are defined such that : 

1 ci 



(25) 



-^4 PU.,, 



r = 



(26) 



where d is a reference length and we introduce coefficient 73 = Ar. The solution 
of system (l25l) is : 

P = P6.+Pl.e-"' + P2.e"^ (27) 
where the boundary conditions at 2; = and yield the values of p^^ and 
P2^- 

(72 + 73)Pi. + (72 - 73)P2, = 7i - 72P6,, 

-e-'^-"(73 - 74)Pi. + e^^"(73 + 74)P2. = -74P6.- 

The liquid density profile is a consequence of Eq. (P7|) when z G [0, /i^;]. 

By taking Eq. (1271) into account in Eq. (1211) and Qb^ip) in linearized form in 

the liquid part of the interlayer, we get immediately 



n(p, 



2 c? 



Pi 



and consequently, 



2cf 



n(PfeJ = — (7i - 72PfeJ(73 + 74) e"*^" + (72 - 73)74Pf,. 
Pi 



X 



(72 + 73)74P6, - (7i - 72PbJ(73 - 74)e^ 



[(72 + 73) (73 + 74)e''"^ + (73 - 74) (72 - 73)e ^^^] 



(28) 



By identification of expressions (1251) and fl2Sl) . we get a relation between 
and Pfj^ and consequently a relation between disjoining pressure n(Pb^) and 
thickness of the liquid film. For the sake of simplicity, we denote finally the 
disjoining pressure by n(/ix) which is a function of at temperature 6. 
Due to the fact that p^^ ^ ~ [16J, the disjoining pressure reduces to the 
simplified expression : 
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2i 

Pi 



(7i - 72P«)(73 + 74)e^'^ + (72 - 73)74P/ 
(72 + 73)74^ - (7i - 72A)(73 - 74)6"''"^ 



[(72 + 73) (73 + 74)e''^^ + (73 - 74) (72 - 73)e"''^1' 



(29) 



Let us notice an important property of a mixture of a van der Waals fluid 
and a perfect gas where the total pressure is the sum of the partial pressures 
of components [45]: at equilibrium, the partial pressure of the perfect gas is 
constant through the liquid-vapor-gas layer -where the perfect gas is dissolved 
in the liquid. The disjoining pressure of the mixture is the same than for a 
single van der Waals fluid and calculations and results are identical to those 
previously obtained. 



6 Numerical calculations for water wetting a vertical plane wall 



Our aim is not to propose an exhaustive study of the disjoining pressure of 
water for all physicochemical conditions associated with different walls but to 
point out examples such that previous results provide new values of maximum 
height for a vertical water film damping a plane wall. 

Calculations are made with Mathematica^^^ . The disjoining pressure 11 and 
the Gibbs free energy G are presented as functions of h^. The hx values must 
be greater than the molecular radius of water corresponding to the smallest 
thickness of the liquid layer. 

The graphs of Tl(hx) are directly issued from Rel. (129!) with different physical 
values obtained in the literature. The graphs of G{hx) are deduced from Rel. 
([3]). As a function of h^, function n(/ix) is not an analytically integrable ex- 
pression; consequently G-graphs are computed by Mathematical but with 
help of a numerical process. 

For a few nanometer range, the film thickness is not exactly hx', at this range we 
must add to the liquid part of the liquid- vapor interface bordering the liquid 
layer (the thickness of which is neglected for films of several nanometers). We 
can estimate this part thickness at 2 ai (half of the thickness of water liquid- 
vapor interface at 20° Celsius [21I3T]). The film thickness is e^^ ~ + 2(Ji. 
The previous results in sections 3-5 remain unchanged by using in place of 
the liquid thickness h. 

When hx = (corresponding to the dry wall), the value of G is the spreading 
coefficient S (see Fig. 2). We must emphasize that point P associated with the 
pancake layer is observed, on the numerical curves, to be closely an inflexion 
point of graph Il{h) corresponding to the strongest stability of films (maximum 
of d'^G/dh'^) To obtain the pancake thickness corresponding to the 
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smallest film thickness, we draw the graphs of Il{hx) and G{hx), when G 
[| ai, i], where £ is a distance of few tens of Amgstrom. 

At = 20° Celsius, we consider successively water wetting a wall in silicon as 
a reference of material and water damping a less wetting wall. 

In c.g.s. units 

The experimental estimates of coefficients are obtained in Refs. [17] and [16] : 

Pi = 0.998, 

Q = 1.478 X 10^ 

cu = 1.4 X 10-^^ 

a I = 2.8 X 10^^ (2.8 Angstrom or 0.28 nanometer), 
mi = 2.99 X 10^23 _ 

From Rel. IQ, we deduce 
A = 1.17 X 10-^ 

72 = 74 = 54.2 . 

From 73 = Ar and Rel. fl26p . we get 

73 = 506, 

d = 2.31 X 10-s. 

We consider two cases : 

a) For silicon (as a reference of wall damped by water), physical character- 
istics are, 

as = 2.7 X 10-^ 
m, = 4.65 X 10-23, 
Psoi = 2.33 . 

No information is available for water-silicon interactions; if we assume that 
cii = cis = 1.4 X 10"^^, we deduce 
7i = 81.2 . 

b) We consider a material such that 7^ = 75 (the material is less damped by 
liquid water). The other values of the coefficients of the material are assumed 
to be the same than in case a). We will see that these values are well adapted 
to our problem. 

Corresponding to Rel. (jl]), the graphs of U^h^c) in cases a) and b) are easy to 
plot following the liquid layer thickness for stable and unstable domains. The 
G-graphs are deduced by numerical integration following the bound h^, the 
limit +00 is replaced by ten thousand molecular diameters of water molecules. 
Due to > |(Ti, it is not possible to obtain numerically the limit point W 
corresponding to the dry wall (Fig. 2). This point is graphically obtained by 
an interpolation associated with the concave part of the C-curve. Point P 
follows from the drawing of the tangent line from W to the G-curve. 
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It is important to point out that reference length d is of the same order than 
ai, as and S and seems a good length order for very thin films. 
In Fig. 3, we present disjoining pressure graphs in the two cases. Real parts 
of disjoining pressure graphs corresponding to dU/dh^ < are plain lines of 
the curves and are associated with thickness liquid layers of several molecules. 
Dashed lines of the curves have no real existence. 

From reporting the pancake thickness of the hx axis for the Il-curve, we de- 
duce its disjoining pressure corresponding value; the maximum of altitude of 
topmost trees is calculated with Eq. (^^. 
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Fig. 3. Graphs represent the disjoining pressure as a function of the liquid layer 
thickness h^ for liquid water at 20° C in contact with different walls (the total 
thickness of the film is Cx ^ h^ + 2o"ij. The unit of x—axis is d = 2.31 x 10^^ 
cm; the unit of y— axis is one atmosphere. Curve (a) corresponds to water wetting a 
silicon wall ('^i = 81.2 cgs); curve (h) corresponds to water wetting a less damped 
wall ('yi = 75 cgs). Points la and lb are the inflexion points of curves (a) and (h) 
corresponding roughly to points P in Fig. 4- 

In Fig. 4, we present graphs of the Gibbs free energy G as a function of hx. The 
limit of the film thickness is associated to the pancake thickness Cp ^ hp + 2 ai 
when the liquid film coexists with the dry wall. The spreading coefficient 
values are associated with point P: at point P, which is numerically close to 
the inflexion point of Il{hx), the tangent goes to point W of the y-axis (where 
G{0) = S'). In the two cases, the total pancake thickness Cp = hp + 2 ai is of 
one nanometer order corresponding to a good thickness value for a high-energy 
surface 

From the graphs, we deduce 5* ~ 64 cgs in case a) and S* ~ 40 cgs in 
case b) corresponding to a less energetic wall. However, crude sap is not pure 
water. Its liquid-vapor surface tension has a lower value than surface tension 
of pure water (72 cgs at 20°C) and it is possible to obtain the same spreading 
coefficients with less energetic surfaces. 
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(a) (b) 




Fig. 4. Graphs represent the film Gibbs free energy per unit area as a function of the 
liquid layer thickness for liquid water at 20° G in contact with two different walls. 
The unit of x— axis is d = 2.31 x 10^^ cm; the unit of y— axis is one cgs unit of 
superficial tension. Graph (a) corresponds to water wetting a silicon wall (ji = 81.2 
cgs); graph (b) corresponds to water wetting a less damped wall (ji = 75 cgs). Point 
W is associated with the surface energy of the dry wall and point P is associated 
with hp corresponding to the pancake layer where the film coexists with the dry wall; 
the smallest film thickness possible is Cp ^ hp + 2a i. 

When is of some hundred meters, Eq. (12^ yields : 

n(a;) ^ piQX. 

The maximum of altitude \Xj^^\ corresponds to the pancake layer. We add 20 
meters to this altitude, corresponding to the ascent of sap due to hydrody- 
namic, capillarity and osmotic pressure. In case h), the material has a lower 
surface energy than silicon and we obtain a film height of 140 meters. 



7 Conclusion 

In this article, we considered a very thin liquid film damping and rising along 
a vertical plane wall as a model of the ascent of sap in xylem tubes. 
At its bottom, the film is bordered on a liquid meniscus corresponding to a 
xylem tube filled with sap up to an altitude of a few ten meters. Above this 
altitude, the xylem tube is embolized such as the liquid water thin film wets 
the wall of the tube. 

The study is static. The motor of the sap motion is induced by the transpira- 
tion across micropores located in tree leaves [3j and is studied by lubrication 
approximation in a model adapted to thin layers pUfTTj . It is natural to fore- 
cast that the diameters of xylem tubes must be the result of a competition 
between evaporation in tubes which reduces the flow of sap and the flux of 
transpiration in micropores inducing the motion strength. 
Computational methods, such as density-functional theory (DFT) and kinetic 
Monte Carlo (KMC) have already had major success in nanoscience. We have 
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obtained the disjoining pressure and Gibbs energy curves for two different 
solid walls. In cases of Lifshitz analysis [18] and van der Waals theory |30j , the 
disjoining pressure behaviors are respectively as 11 ~ and 11 ~ exp{—h). 
These two behaviors seem unable to study a film with a thickness between 
one and three nanometers. 

It is wondering to observe that the density-functional theory expressed by a 
rough model correcting van der Waals' with a surface density-functional at 
the walls enables to obtain a good order of magnitude of the ascent of sap. 
This result is obtained without any complex weighted density-functional and 
without taking into account the quantum effects corresponding to less than an 
Amgstrom length scale. The surface density-functional at the wall takes into 
account a power-law behavior associated with a balance between attractive 
and repulsive forces, the square-gradient functional schematizing the liquid- 
vapor interface effects. This biophysical observation seems to prove that this 
kind of functional can be a good tool to study models of liquids in contact 
with solids at a small nanoscale range. 

Due to a remark by James R. Henderson [18] , it is interesting to note that if we 
switch the micro-tube surfaces to wedge geometry as in [39] or to corrugated 
surface then, it is much easier to obtain the complete wetting requirement. 
Thus, plants can avoid having very high energy surfaces, but still be internally 
wet, if they pass liquids through wedge shaped corrugated pores. The wedge 
does not have to be perfect on the nanometer scale to significantly enhance 
the amount of liquid that would be passed at modest pressures corresponding 
to nano-sized planar films. It is bound to improve on the calculation because 
it enhances the surface to volume ratio. 

In such a case, we remark that the wall boundary can always be considered 
as a plane surface with an average surface energy as in Wenzel's formula |50j . 

Acknowledgment: The author dedicates the paper to the memory of Professor 
Pierre Casal, his Master and friend. 
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